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Abstract 

We study J/ip production in pp collisions at RHIC within the PHENIX detector acceptance 
range using the color singlet and color octet mechanism which are based on pQCD and NRQCD. 
Here we show that the color octet mechanism reproduces the RHIC data for J /ip production in pp 
collisions with respect to the px distribution, the rapidity distribution and the total cross section 
at ^fs = 200 GeV. The color singlet mechanism leads to a relatively small contribution to the total 
cross section when compared to the octet contribution. 
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Understanding the J /if) production mechanism in high energy hadronic collisions is an 
important topic in QCD and in collider physics. Two prominent mechanisms for J /if) pro- 
duction at high energy colliders are l)the color singlet mechanism and 2)the color octet 
mechanism both of which rely on pQCD and NRQCD. At Tevatron energies, the color 
singlet mechanism was found to give too small a yield and the color octet mechanism for J/ if) 
production was then introduced 0,0]. It was shown that the color octet mechanism, with 
both parton fusion and parton fragmentation processes included, reproduced the Tevatron 
data |6||. The relativisitic heavy ion collider (RHIC) at Brookhaven is a unique facility which 
can collide unpolarized and polarized protons at y/s =200 and 500 GeV as well as two heavy 
nuclei such as gold on gold. It has also the capability of exploring light-heavy interactions 
such as the collisions pAu and dAu at \fs = 200 GeV. The unpolarized pp collision data 
serve as a baseline for understanding collective phenomena such as the production of a quark 
gluon plasma that is expected in the heavy ion collisions as well as providing the necessary 
information needed to extract interference terms generated when considering polarized pp 
collisions. One of the key things to be measured in Au Au collisions as well as in polarized 
pp collision is J /if) production. In relativistic heavy ion collisions, J/ if) suppression has been 
posited to be one of the major signatures of the production of the quark-gluon plasma j^. 
In case of polarized pp collisions the measurement of J /if) production is an important in- 
gredient in extracting the polarized gluon distribution function for the proton directly. For 
that reason, understanding the dominant mechanism of J /if) production is important. Thus 
it is crucial to first determine the relative importance of color singlet and color octet con- 
tributions in the existing unpolarized proton-proton scattering data at RHIC energies. In 
this paper we will study J / if) production in pp collisions at RHIC by using the color singlet 
and color octet mechanism within pQCD and NRQCD. We will study the pt distribution, 
rapidity distribution and total cross section of the J /if) production at yfs = 200 GeV pp 
collisions and compare those with the recent measurements by the PHENIX collaboration 
at RHIC J 

In the color singlet mechanism, quarkonium is formed as a non relativistic bound state of 
a heavy quark- ant iquark pair via static gluon exchange. In this mechanism it is assumed that 
the QQ pair is produced in the color singlet state at the production point with appropriate 
spin (S) and orbital angular momentum quantum number (L) and then evolves into a bound 
state 2S+1 Lj with total angular momentum quantum number (J). The relative momentum 
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of the QQ pair inside the quarkonium is assumed to be small compared to the mass m of 
the heavy quark so that the Q and Q will not fly apart to form heavy mesons. The non- 
relativistic non-perturbative wave functions and its derivatives appearing in the color singlet 
calculation are either determined from the potential model or taken from experiment. 

The leading term in the Pt distribution for the cross section for J/ip production first 
occurs at order oP s . The pt distribution for heavy quarkonium production in the color 
singlet mechanism is given by: 

^-(ab -> j/?/>, X j + x) = Y,J d yJ dx * x * 

fa/A{Xa,Q 2 ) x'j b/B (x' b ,Q 2 ) ^ 

| (ab - C^Lj)c) , (1) 

where 

1 x b ^SM T e-y - M 2 
6 x b ^-M T ev ■ 1 ] 

In the above equation a, b, c are light quarks and gluons. f(x, Q 2 ) is the parton distribution 

function at longitudinal momentum fraction x and at factorization scale Q. M is the mass 



of the bound state quarkonium and M T = \JPt + M 2 . The partonic level differential cross 
section ^(&b — * ( 2S+1 Lj)c) contain the non-relativistic wave function |-R[0)| 2 (for direct 
J/ip process) and its derivatives |_R'(0)| 2 (for xj processes) at the origin 



For the non 



relativistic wave functions at the origin we take the Buchmuller-Tye wave function 91. The 

: |i?(0)| 2 =0.81 GeV 3 . For the derivative of the radial wave function at 
origin we use |lo|: ^^2^=15 MeV. 

In the color octet mechanism the relativistic effects are taken into account which are 
neglected in the color singlet case. In the color octet mechanism using an effective field 
theory called non-relativistic QCD (NRQCD), the dynamical gluon enters into Fock state 
decompositions of the quarkonium states. In NRQCD the expansion is carried out in terms 
of the relative velocity v (v 2 ~ 0.23 for CC system and 0.1 for BB system) of the QQ bound 
state. The NRQCD Lagrangian density is given by 

^NRQCD flight ^heavy ^correction- (3) 

In the above equation light refers to light quarks and gluons, heavy refers to lowest order 
non-relativistic heavy quarks part and correction is higher order corrections in heavy quark 



sector. The explicit expressions for the three terms can be found in Q]. In NRQCD the 
dynamical gluons enter into the Fock state decompositions of different physical states. The 
wave function of an S-wave orthoquarkonia state \ipQ > is expanded as follows: 

> = 0(1)\QQ[ 3 S[ 1} ] > +0(v)\QQ[ 3 pf]g > 
+0(v 2 )\QQ[ 3 S[ 1 %g > +0(v 2 )\QQ[ l Si% > 

+0(v 2 )\QQ[ 3 D { j%g>+ (4) 

and the wave functions of a P-wave orthoquarkonium state \xqj > has a similar expansion: 

\x Q j >= 0(1)\QQ[ 3 P^} > +0(v)\QQ[ 3 s[% > 

+0(v 2 ).... (5) 

In the above equation (1,8) refers to singlet and octet state of the QQ pair. After a QQ 
is formed in its color octet state it may absorb a soft gluon and transform into \xqj > via 
eq. (0) and then become a J/ip by photon decay. The QQ pair in the color octet state 
can also emit two long wavelength gluons and become a J/ip via eq. (J1J) and so on. The 
non-perturbative matrix elements can be fitted from other experiments or can be determined 
from lattice calculations. 

In the color octet mechanism the distribution differential cross section for J/ijj pro- 
duction is given by: 

da 



dp> 



T 



(AB -> 4>q(xj) +X) = j d V J ' dx " x * 



x a - 



fa/A{x a ,Q ) x' b f b/B (x' b ,Q ) X 



§(ab^CC[ 2S ^Lf]c^M Xj )), (6) 
at 



where the partonic level differential cross section is given by 
% {ab^CC[^Lf]c^^ Q { X j)) ' 



£ \A(ab - CC^Lf ]c) short | 2 < 0\Ot Q{xj \ 2S+1 Lj)\0 > (7) 
In our calculation all the partonic level matrix elements squared 

£ \A(ab^CC[ 2S+1 Lf]c) shoTt \ 2 
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and the non-perturbative matrix elements 

<0\oi Q{xj \ 2S+1 Lj)\0> 

are taken from 

For the rapidity distribution and total cross section of J /if) production the leading con- 
tribution is proportional to a 2 and we will only consider that contribution here. In leading 
order, the parton fusion processes contribute to the total cross section for J /if) (xj) produc- 
tion in pp collisions as follows: 

a PP^J/4>(xj) = Sa fe J dXa J dx b f a/A (x a ,Q 2 )f b/B (x b ,Q 2 ) 

a ab-.JWxj) (s)5(x a x b - Am 2 /s) (8) 

where a ab ^ J/ ^ XJ (s) is the partonic level cross section at leading order which are derived in 



[5, In the above equation m is the mass of the charm quark. At the leading order not 
all the processes contribute to the total cross section in the singlet and octet channel. In 
the color singlet channel the J /if) production cross section at a 2 order is given by: 

<tir J '%) = oT*\s)BR X0 , + a^(s)BR X2 . (9) 

Here BR X refers to the branching ratio of a x to decay into a J /if). Similarly in the octet 
channel at the leading order the J / if) production cross section is given by: 

+air xo (s)BR X0 + a^ X2 (s)BR X2 . (10) 

In the above equations the direct J /if) production and the decay of Xc with their correc- 
sponding branching ratios to J / if) are included. In the leading order we use the same values 
of the non-perturbative matrix elements used in the singlet and octet channels as in 

We will use the GRV980 and MRST99 parton distribution functions inside the 
proton. The factorization scales are set to be Q = 2m for the total cross section and 



rapidity distribution plots, and Q = \jv\ + M 2 for the plot of the px distribution. 

In Fig.l we present the pt distribution of the J /if) production cross section separately 
for the color singlet and color octet mechanism including the effects of cuts due to the 
°HENIX detector acceptance range. We have considered all the parton fusion processes 
5, 13] in our calculation. In these plots we have used the GRV98 and MRST99 parton 



distribution functions. As the non-perturbative NRQCD matrix elements has uncertainties 
as given in we present two plots for each of the PDF's which correspond to the upper 
and lower limit on the non-perturbative matrix elements. The solid line corresponds to 
the results obtained by usingthe GRV98 PDF's with lowest values of the non-perturbative 
matrix elements as given in |5J]. The upper dashed line corresponds to the results obtained 
by using the GRV98 PDF's with highest values of the non-perturbative matrix elements. 
The upper dot-dashed line corresponds to the results obtained by using the MRST99 PDF's 
with highest values of the non-perturbative matrix elements. The lower dot-dashed line 
corresponds to the results obtained by using the MRST99 PDF's with lowest values of the 
non-perturbative matrix elements. The lower dashed line corresponds to the color singlet 
contribution. We have compared our calculation with the recent PHENIX data (run3). The 

n 

daggers are recent PHENIX data taken from (SJ). It can be seen from the figure that the color 
singlet contriubtion is more than an order of magnitude smaller than the octet contribution. 
On the other hand, the color octet mechanism with only parton fusion processes reproduces 
the recent PHENIX data very well. Thus for p? < 5GeV we do not find it necessary 
to include contributions from the parton fragmentation processes. We have not compared 
our calculation with the data below pt equals to 2 GeV as pQCD calculations are not 
valid at such low transverse momenta. In the near future PHENIX will have sensitivity up 
to pt = lOGeV. This new data will allow us to decide whether fragmentation processes 
contributions within color octet mechanism are important at higher p? and whether their 
inclusion will allow us to understand the data. 

In Fig. 2 we present the rapidity distribution of the J/ip production cross section in pp 
collisions at RHIC at = 200 GeV. We have taken the charm quark mass equal to 1.45 GeV 
in our calculation. The solid line is the results obtained by using the GRV98 PDF's with 
lowest values of the non-perturbative matrix elements. The upper dashed line corresponds to 
the results obtained by using the GRV98 PDF's with highest values of the non-perturbative 
matrix elements. The upper dot-dashed line corresponds to the results obtained by using the 
MRST99 PDF's with highest values of the non-perturbative matrix elements. The lower dot- 
dashed line corresponds to the results obtained by using the MRST99 PDF's with lowest 

n 

values of the non-perturbative matrix elements. The PHENIX experimental data |8J are 
shown in the figure as daggers. The rapidity range covered by PHENIX at RHIC is from 
—3 < y < 3. The color singlet contribution is also shown as the lower dashed line (here 
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we only show the results using the GRV98 PDF's) which is very small when compared with 
the color octet contribution. It can be seen that just including the color octet contribution 
explains the run3 PHENIX data. 

Finally, we present the total cross section for J/ip production in pp collisions at y/s = 
200 GeV at RHIC and compare it with the run3 PHENIX data. In the color octet case we 
find Baji^— 159 ± 6 (nb) which is in very good agreement with the PHENIX collaoration 
experimental value which is reported as 159 (nb) jsj]. The ± 6 (nb) is due to the errors 
in the non-perturbative matrix elements as given in [^|. In the above total cross section 
computation we have used charm quark mass equal to 1.45 GeV and the MRST99 parton 
distribution function inside the proton. The agreement of color octet mechanism predictions 
with the PHENIX experimental data is remarkable in all aspects, i.e. in distribution, in 
rapidity distribution and in total cross section measurements at y/s = 200 GeV pp collisions. 

Before summarizing, we briefly discuss the uncertainties due to PDF's, non-perturbative 
matrix elements etc. As can be seen from Fig. 1 and Fig. 2 that both GRV98 and MRST99 
PDF's describe the data very well. We have also checked that the recent CTEQ6M PDF's 
also explain the data very well. The uncertainty in the results due to the errors in the 
non-perturbative matrix elements are given in Fig. 1 and 2 and also in the total cross 
section. 

In summary, in this letter, we have studied J/ip production in pp collisions at RHIC at 
y/s = 200 GeV within the PHENIX detector acceptance range by using the color singlet and 
the color octet mechanism within pQCD and NRQCD and have compared them with the 
recent PHENIX experimental data. We have demonstrated that the color octet mechanism 
using only parton fusion processes is able to reproduce the RHIC data on the distribution, 
rapidity distribution and total cross section of J ftp production. This is done without any 
normalization factor modifications and using matrix elements extracted from the Tevatron 
and other experiments. The color singlet mechanism gives negligible contributions when 
compared to the octet contributions at this energy. Our results confirm the applicability of 
the color octet mechanism found at the Tevatron. 

It is important to identify the correct mechanism for J/%p production in pp collisions at 
RHIC if we want to study J/ ip suppression as a signature of quark-gluon plasma detection in 
AuAu collisons at y/s = 200 GeV 7 1 . It is also important if we hope to extract the polarized 
gluon distribution function from the data on polarized pp collisions at RHIC. The latter will 
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be important for isolating interference terms which can be used to explore various models 
which predict non standard model TeV physics. 
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FIG. 1: pt distribution of J /if) production cross section in pp collisions at yfs = 200 GeV at RHIC 
in the PHENIX detector acceptance range. The solid and upper dashed lines and upper/lower dot- 
dashed lines correspond to color octet mechanism predictions with GRV98 and MRST99 PDF's. 
The lower dashed line is the color singlet contribution. The daggers are run3 PHENIX experimental 
data. 
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FIG. 2: Rapidity distribution of J /if) production cross section in pp collisions at y/s = 200 
GeV at RHIC in the PHENIX detector acceptance range. The solid and upper dashed lines and 
upper/lower dot-dashed lines correspond to color octet mechanism predictions with GRV98 and 
MRST99 PDF's. The lower dashed line is the color singlet contribution. The daggers are run3 
PHENIX experimental data. 
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